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We develop a graphical calculus for completely positive maps and in doing so cast the theory of 
open quantum systems into the language of tensor networks. We tailor the theory of tensor networks 
to pictographically represent the Liouville-superoperator, Choi-matrix, process-matrix, Kraus, and 
system-environment representations for the evolution of open-system states, to expose how these 
representations interrelate, and to concisely transform between them. Several of these transfor- 
mations have succinct depictions as wire bending dualities in our graphical calculus — reshuffling, 
vectorization, and the Choi-Jamiolkowski isomorphism. The reshuffling duality between the Choi- 
matrix and superoperator is bi-directional, while the vectorization and Choi-Jamiolkowski dualities, 
from the Kraus and system-environment representations to the superoperator and Choi-matrix re- 
spectively, are single directional due to the non-uniqueness of the Kraus and system-environment 
representations. The remaining transformations are not wire bending duality transformations due 
to the nonlinearity of the associated operator decompositions. Having new tools to investigate old 
problems can often lead to surprising new results, and the graphical calculus presented in this paper 
should lead to a better understanding of the interrelation between CP-maps and quantum theory. 



I. INTRODUCTION 

A complete description of the evolution of quantum 
systems is an important tool in quantum information pro- 
cessing (QIP). In contrast to closed quantum systems, in 
open quantum systems the evolution need no longer be 
unitary. In general the evolution of an open quantum 
system is called a quantum operation or quantum chan- 
nel which, for a discrete time interval, is described math- 
ematically by a completely positive map (CP-map) [T]. 

In the context of QIP, a quantum channel is a com- 
pletely positive linear map acting on the density opera- 
tors which describe a physical systems state. A map £ 
is positive if and only if it preserves the positivity of an 
operators spectrum, and completely positive if and only if 
it further satisfies the condition that the composite map 
T® £ is positive, where I is the identity map on a space 
of density operators of at least equal dimension as the 
space on which £ acts. We will consider the case when £ 
is also trace preserving, which is to say Tr[£(p)] = Tr[p] 
for all density operators p. Since quantum systems are 
described by density operators, positive operators p with 
unit trace, the requirements that £ be complete-positive 
and trace-preserving ensures that our output state will 
always be a valid density operator. Such maps are called 
completely positive trace preserving maps (CPTP-maps). 

There are numerous representations for completely- 
positive maps, and while these representations are at 
least in principle equivalent, transformations between 
representations is often cumbersome and tedious. Thus 
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a key result in the present paper is to introduce and de- 
velop graphical calculus methods that facilitate an intu- 
itive unification and interoperability between these rep- 
resentations, the outcome of which is depicted in Fig. [T] 

The graphical calculus casts the theory of open quan- 
tum systems into the framework of tensor networks, 
which comes equipped with a graphical means to rep- 
resent and reason about the contraction of sequences of 
tensors [2]. Graphical calculi have been used to great 
benefit in several areas of modern physics with the most 
prolific example being the use of Feynman diagrams to 
calculate scattering amplitudes in quantum field theo- 
ries [5]. In the case of tensor networks, their use dates 
back to earlier work by Penrose who's graphical nota- 
tion is a useful starting point [2]- Recent work has been 
done on uniting the theory of tensor network states with 
modern algebra and similar graphical tools which have 
recently received interest in the quantum information, 
foundations and condensed matter communities [4j [5]. 
They have been used as computational tools for simulat- 
ing many-body quantum systems efficiently [4, 6-5], as a 
tool for manipulating tensor networks and to generalize 
quantum circuits [5j, and for exploring the foundations 
of quantum mechanics [5HT2]. 

Although it is straightforward to translate equations 
into so called tensor string diagrams, a missing piece has 
been a graphical calculus for open systems theory which 
provided new results, and hence enhanced the potential 
for diagramatic reasoning. Several proposed graphical 
languages have been based on category theory or a gen- 
eralization of quantum circuits (SHS 1131 114] . Our work 
differs from the current approaches. We develop a cal- 
culus and tailor the tensor string diagrams of Penrose to 
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FIG. 1. The main mathematical representations for completely positive maps and how one may transform between them. Solid 
arrows represent linear operations which we prove can be done by "wire bending" transformations in our graphical calculus. 
Dashed arrows represent non-linear transformations. Reshuffling and Stinespring dilation are bijective transformations, vec- 
torization and the Jamiolkowski isomorphism are surjective transformations, and the spectral decomposition is an injective 
transformation. 



unify several mathematical representations used in open 
quantum systems and to transform freely between them. 
In accomplishing this, we express our graphical tensor 
calculus in the Dirac notation familiar in QIP instead of 
the abstract index notation used by Penrose. We have 
also introduce a graphical colour summation convention 
based on colouring of tensors in our tensor networks to 
represent summation over indices. This appears to ex- 
tend the capabilities of diagrammatic reasoning, with po- 
tential applications to our techniques, and semantics even 
outside of quantum physics. 

The present paper has three sections, in Section |IT| we 
introduce the elements of our graphical calculus with a 
particular emphasis on bipartite systems and vectoriza- 
tion of linear operators that are important in the subse- 
quent sections. In Section |III| we introduce the math- 



ematical representations of CP-maps listed in Fig. [T] 
and their corresponding graphical representation. In Sec- 
tion [IV] we describe how one may transform between any 
of the representations of CP-maps as shown by the ar- 
rows in Fig.[l] In this figure the type of arrow connecting 
two boxes depicts the type of transformation between the 
corresponding CP map representations. The solid arrows 
depict linear transformations, while the transformations 
depicted by dashed arrows are non-linear. Here the util- 
ity of our graphical calculus manifests as all the linear 
transformations correspond to wire-bending operations 
in our tensor diagrams, yet the corresponding algebraic 
expressions involve complicated manipulations. In ad- 
dition, the bi-directional arrows depict invertible (bijec- 
tive) transformations, while the single directional arrows 



depict non-invertible transformations, where either injec- 
tivity (solid arrows), or surjectivity (dashed arrow) fails. 

In the present paper we will use the notation that 
X,y, Z are finite-dimensional complex Hilbert spaces, 
£(X, y) is the space of bounded linear operators A : X — > 
y (with C(X) = C(X, X)), T(X,y) is the space of op- 
erator maps £ : £{X) C(y) (with T{X) = T(X,X)), 
and C(X,y) is the space of operator maps £ which are 
CPTP 



II. TENSOR NETWORKS 

Tensors can bethought of as indexed multi-dimensional 
arrays of complex numbers with respect to a fixed stan- 
dard basis. The number of indices is called the order 
of a tensor, and the concurrent evaluation of all indices 
returns a complex number. For example, consider the 
Hilbert space X = C d , where as is typical in QIP, we 
choose our standard basis to be the computational ba- 
sis {\i) : i — 0, ...,d — 1}. Then in Dirac notation a 
vector \v) £ A" is a lst-order tensor which can be ex- 
pressed in terms of its tenor components Vi := (i\v) with 
respect to the standard basis as \v) = X2j=o Sim- 
ilarly one can represent linear operators on this Hilbert 
space, A £ £(X), as 2nd-order tensors with components 
A ij :={i\A\j) aa A = ^-i A ij \i){j\. 

Hence, in Dirac notation the number of indices of a ten- 
sors components are what we refer to as the order of the 
tensor. Vectors \v) £ X refer to tensors which only have 
ket basis elements, vectors in the dual vector space 
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(u\ e X^ refer to those with only bra's "{if, and linear 
operators on A € /^(A") refer to tensors with a mixture 
of kets and bras in their component decomposition. 

The idea of representing states, operators and process 
(etc.) diagrammatically dates back to several works by 
Penrose and is often referred to as Penrose graphical no- 
tation or string diagrams. We adopt Penrose's notation 
of representing states (vectors) and effects (dual-vectors) 
as triangles, linear operators as boxes, and scalars as di- 
amonds, as illustrated in Fig. [2j Here each index corre- 
sponds to an open wire on the diagram and so we may 
define higher order tensors with increasingly more wires. 
The number of wires is then the order of the tensor, with 
each wire acting on a separate vector space Xj. 
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for fixed dimension d, though one may generalize most of 
what follows to situations where the dimensions of each 
wire are not equal. Note that we represent scalars as ten- 
sors with no open wires, this could either be a contracted 
tensor A = (v\u), or a multiplicative factor A acting on 
the tensor A as XA. 
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(a) Vector transposition: 
\v) T = (v\ 



(b) Dual-vector 
transposition: (?;| T = \v) 



(c) Linear operator transposition 
FIG. 3. Graphical depiction of transposition by wire bending 
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FIG. 2. Elementary tensors. 

We also insist that the orientation of these wires, rather 
than the number of wires, specifies whether they repre- 
sent multi-partite vectors, dual-vectors, or linear opera- 
tors. We have a freedom in choosing our orientation for 
the tensors, top-to-bottom, bottom-to-top, left-to-right 
or right-to- left. In this paper we will choose the right-to- 
left convention (the opposite of most orthodox quantum 
circuits) so that the graphical representation will most 
closely match the underlying equations. Thus we use the 
terms vector, dual-vector and linear operator to refer to 
tensors of any order, not just lst-order and 2nd-order, 
based on the orientation of their wires as follows: 

1. Vectors \v) G ®i=i are tensors with n wires 
oriented to the left. 

2. Vectors in the dual space (v\ £ <S)™ = i are tensors 
with n wires are oriented to the right. 

3. Linear operators A : (£>™ =1 Xi — > <S) J= i are ten- 
sors which have n wires going to the right and m 
wires to the left. 

4. Tensors with no open wires are Scalars A E C. 

The graphical depictions of the these tensors are also il- 
lustrated in Fig. [2] In the present paper, we will generally 
be interested in the case where each wire indexes Xa = C d 



The mathematical rules of tensor network theory as- 
sert that the wires of tensors may be manipulated, with 
each manipulation corresponding to a specific operation 
with equational meaning. We now introduce some tools 
which we have tailored for manipulations common in 
open quantum systems. Transposition of a 2nd-order 
tensor is represented a bending of a tensors wires as il- 
lustrated in Fig. [3] Complex conjugation of a tensor's 
coefficients however is depicted by a bar over the tensor 
label in the diagram as illustrated in Fig. [4] Hence we 
may represent the transformation of a vector to its dual 
vector, or the hermitian conjugation of a linear operator 
as the combination of these two operations as illustrated 
in Fig. [5] We stress that under this convention a vector 
\v) = ^2.iV l \i) and its hermitian conjugate dual-vector 
(v\ = J^i are represented as shown in Fig 
respectively. 



2a 



and 2b 
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(a) Complex 
conjugation of 
vector \v) 



(b) Complex 
conjugation of 
dual vector (v\ 



(c) Complex 
conjugation of 
operator A 



FIG. 4. Graphical depiction complex conjugation. 

Tensor contraction is represented by joining the wires 
corresponding to the indices to be contracted. In the case 
of matrix multiplication A ■ B is represented by connect- 
ing the corresponding wires of the tensors representing 
the matrices as shown in Fig. 



Ga 



To form multi-partite 
tensors we denote the tensor product of two tensors A®B 
by the vertical juxtaposition of their tensor networks as 
shown in Fig. 6b The trace of an operator by depicted 
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(a) Hermitian conjugation (b) Hermitian conjugation 
of a vector : = (v\ of a dual- vector : (v\^ = |u) 




(c) Hermitian Conjugation of an operator A 
FIG. 5. Graphical depiction of hermitian conjugation. 




FIG. 8. Illustration of our colour summation convention using 
the spectral decomposition of a normal operator. 

■■= |X3— 

FIG. 9. Graphical resolution of the identity using the colour 
summation convention. 



by connecting the corresponding left and right wires of a 
linear operator as illustrated in Fig. [7] 
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(a) Matrix multiplication 
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(b) Tensor product 

FIG. 6. Vertical (b) and horizontal (a) Composition of ten- 
sors. 
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FIG. 7. Tensor network for the trace of a linear operator A. 



The unnormalized maximally entangled Bell-state 

d 

\&+)=Y^\£)®\i)eX®x, (l) 
«=o 



is represented graphically as the curve shown in Fig. |10a] 
Similarly the Bell-effect ($ + | is represented as shown in 
Fig. [10b} 



= 1> 



(a) Bell-state 
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(b) Dual Bell-state or Bell effect 

FIG. 10. The graphical depiction of the unnormalized Bell- 
state |$ + ) and its dual ($ + |. 



We represent summation in our variant of graphical 
calculus adapted to open systems by introducing shad- 
ing or colouring of the tensors being summed over (we 
call this the colour summation convention). Tensors cor- 
responding to the same summation index will be shaded 
the same colour, and we use different colours for different 
summation indexes. For example, consider the spectral 
decomposition of a normal operator A with eigenvalues 
and eigenvectors A* and |af) respectively. The graphical 
spectral-decomposition A = J2i^i\ a i) using our graphi- 
cal summation convention is illustrated in Fig. [5J 

In our colour summation convention, we will represent 
the sum over the standard basis as a shaded vector (or 
dual-vector) tensor with an empty label. This is illus- 
trated in the graphical resolution of the identity appear- 
ing in Fig [9] 



As will be shown in Section |"HB| our choice of graph- 
ical notation for \$ + ) is due to its equivalence to the 
vectorization of the identity operator. 

Using the graphical definition for |$ + ) we can compose 
the bell-state and its dual to form an identity element 
[2] . This is known as the snake equation and is shown in 

Fig. OH 

The snake equations have several uses and provide an 
equivalence class of diagrams. Anytime we have a curved 
wire with two bends we can "pull the wire" to straighten 
it out into an identity. Anytime we bend a wire, trans- 
forming between say a bra and a ket, we can bend the 
wire to transform back again. 

By combining the snake-equation with the wire- 
bending operation for transposition from Fig.|3cJ we find 
that "sliding" a linear operator around a bell-state is also 
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FIG. 11. The snake or zig-zag equations [2]. 



equivalent to transposition of the operator. This is illus- 
trated in Fig. [l2j Note that due to the orientation of the 
wires this graphical representation of the operator A is 
actually a vector. This is called the vectorization of a 
matrix and we discuss this in more detail in Section [II Bl 




FIG. 12. Graphical depiction of transposition by "sliding" 
around a bell-state. 



A. Bipartite Matrix Operations 

Bipartite matrices are used in several representations 
of CP-maps, and manipulations of these matrices will 
be important in the following discussion. Consider two 
complex Hilbert spaces X , and y with dimensions d x and 
d y respectively. The bipartite matrices we are interested 
in are then d 2 x x dy matrices M G C{X ® y) which we can 
represent as 4th-order tensors with tensor components 

M m ^ nv := (m,fi\M\n,v) (5) 

where m, n G {0, ...,d x — 1}, (X, v € {0, ...,d y — 1} and 
1 71, v) :— \n) ® \v) € X Cgi y is the tensor product of the 
standard bases for X and y. 

We can also express the matrix M as a 2nd-order ten- 
sor in terms of the standard basis {\a) : a = 0, . . . , D — 1} 
for X eg) y where D — d x d y . In this case M has tensor 
components 



Another important operation is the graphical SWAP 
which exchanges the position of two Hilbert spaces in 
a composite system. Let X and y be complex Hilbert 
spaces of dimensions d\ and d 2 respectively, then the 
SWAP operation is the map 



swap : x <g> y -> y <gi x 

SWAP : |x) <g> |y> h+ |y) \x), 



(2) 
(3) 



for all \x) e X, \y) € y. 

Given any two orthonormal basis {\xi) : i = 0, . . . , d\ — 
1} and {\yj) : j = 0, . . . , d% — 1} for X and y respectively, 
we can give an explicit construction for the SWAP oper- 
ation as 



SWAP = ^2 £ \y J ){x l \®\x l ){y 3 



(4) 



i 1= j 2 =0 



The SWAP operation is represented graphically by two 
crossing wires as shown in Fig. 13 The basis decomposi- 
tion in Eqn Q is then an application of the resolution of 
the identity from Fig. [9] to each wire. In Section |HB| we 
will see that the SWAP operation is the natural trans- 
formation between the row-stacking and column-stacking 
vectorization conventions. 




x 



FIG. 13. The SWAP operation in our graphical calculus (left) 
and it's decomposition in terms of the standard basis using 
our colour summation convention (right). 



M, 



a/3 



(a\M\ 



(6) 



This is represented graphically in Fig. [14} We can specify 
the equivalence between the tensor components M a p and 
M m(t)n „ by making the assignment 

a = dyin + n 
f3 = d y n + v, 

where d y is the dimension of the Hilbert space y. 



(7) 
(8) 



M, 



m/i/iv 




(a) 4th-order tensor components of a 
bipartite matrix M 




(b) 2nd-order tensor components of a bipartite matrix M 

FIG. 14. Tensor components of a Bipartite matrix M G C(X® 
X). 



The bipartite matrix operations which are the most 
relevant for open quantum systems (see Fig. [IJ are the 
partial trace over X (Tr^) (and Try over y), transposi- 
tion (T), bipartite- SWAP (S), col-reshuffling (R c ), and 
row-reshuffling (R r ). The corresponding graphical ma- 
nipulations are in Fig. [15] and in terms of the tensor com- 
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ponents of M, these operations are respectively given by: one of two mappings 



Tr* : C{X ® y) — ► c(y) 



7nfi,rnv 



Try : £(# <g> y) — > £(#) 
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s ■. c{x®y) — ► c{y®x) 



Mr, 



M„ 



r c : c(x ®y)^ c(x ®x,y®y) 



m„ 



M„ 



R r : c(x ® y) — > c(y ®y,x® x) 



M. 



7rtfj,,nv 



Note that we will generally use reshuffling R to refer to 
col- reshuffling R c . 




(a) Partial Trace (b) Partial Trace 
Tr^(M) Try(M) 



(c) Transpose 
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(d) Bipartite-Swap (e) Row-reshuffle (f) Col-reshuffle 



M 1 



M 1 



FIG. 15. Tensor networks for common Bipartite matrix ma- 
nipulations. 



B. Vectorization of Matrices 

We now recall the concept of vectorization which is a 
reshaping operation, transforming a (mxn)-matrix into a 
(1 x mn)-vector |15j . This is necessary for the description 
of open quantum systems in the superoperator formal- 
ism, which we will consider in Section [lII C| Vectorization 
can be done with using one of two conventions: column- 
stacking (col-vec) or row-stacking (row-vec). Consider 
two complex Hilbert spaces X = C m ,y = C™, and linear 
operators A £ £(X,y) from X to y. Vectorization is 



row-vec: C(X, y) -> y ® X 
col-vec: C(X,y) -J- X ® y 



A^\A)) r (9) 
A^\A)) C (10) 



where the operation col(row)-vec when applied to a ma- 
trix, outputs a vector with the columns (rows) of the 
matrix stacked on top of each other. Graphical represen- 
tations for the row-vec and col-vec operations are found 
from bending a wire to the left either clockwise or coun- 
terclockwise respectively. This is illustrated in Fig. |16| 



— A 



(a) Row-vec 




(b) Col-vec 

FIG. 16. Tensor networks for Vectorization in the (a) row 
and (b) column conventions are found from a wire bending 
duality. 



Vectorized matrices in the col-vec and row-vec conven- 
tions are naturally equivalent under wire exchange (the 
SWAP operation) as illustrated in Fig. 17 




FIG. 17. Col-vec and row-vec vectorized matrices are natu- 
rally equivalent under the SWAP operation. 

We may also define a vectorization operation with 
respect to an arbitrary basis for C(X, y). Let X = 
C4 , y — C dy , and D = d x d y . Vectorization with respect 
to an orthonormal operator basis {a a : a = 0, D — 1} 
for £(X,y) is given by 



cr-vec: C(X) -> X ® X : A h> \A)) a . 



(11) 



This operation extracts the coefficients of the basis ele- 
ments returning the vector 



D-l 



\A)) a :=X>[atA]|a> 



(12) 



a=0 



where {\a) : a = 0, ...,D — 1} is the standard basis for 
X (g) y. This is depicted in our graphical calculus in 
Fig. [18} 

To distinguish between these different conventions use 
the notation \A)) X to denote the vectorization of a matrix 
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FIG. 18. Tensor network for a- vectorization ( 12 1 with respect 
to an orthonormal operator basis {a a }. 



A, were the subscript x — c,r,o labels which convention 
we use; either c for col-vec, r for row-vec, or a for an 
arbitrary basis. 

For the case X = y = C d , we can define row-vec and 



col-vec in terms terms of Eqn ( 12 ) by taking our basis 
to be the elementary matrix basis {Ei r j — : i,j = 

0, d 2 — 1} and making the assignment a = di + j and 
a = i + dj respectively. Hence we have 



d-l 

\A))r ■■= E ^ I*) ® \fi 
i,j=0 
d-l 

\A)) C := £ ^ li) ® \i) 

i,j=0 



(13) 



(14) 



for row-vec ( [L3| and col-vec (14 1. From the dehnition of 
the Bell state |$ + ) and summing over i and j in Equa- 
tions ( 13 1 and ( 14 ) respectively one can rewrite this as 



\A)) r = (4g 1)|$ H 
|A)) C = (1®A)\& 



(15) 
(16) 

which are the equational versions of Fig.s |16a| and |16b| 
respectively. 

When working in the superoperator formalism for open 
quantum systems, it is sometimes convenient to trans- 
form between vectorization conventions in different bases. 
Given two orthonormal operator bases {<J a } and {w Q } for 
C(X,y), the basis transformation operator 

: X®y^ X®y : \A)) a ^ \A)) U (17) 

transforms vectorized operators in the cr-vec convention 



to the w-vec convention as illustrated in Fig. 19 





FIG. 19. Tensor network for the change of vectorization basis 
for a vectorized operator A from the cr-vec convention to the 
oj-vec convention. 

The basis transformation operator T^^^ is given by 
the equivalent expressions 



=J2\a)((uj a \ <T 

a 



(18) 
(19) 



and the corresponding tensor networks are given in 
Fig. [20] The proof's of Eqn ( |18|19[ ) are found in Ap- 



pendix IB 




FIG. 20. Tensor networks for the change of vectorization basis 
operator T^^u from the cr-vec to w-vec conventions, where 
{<r„}, {u} a } are orthonormal bases for C(X,y). 

For the remainder of this paper we will use the col-vec 
convention by default, and drop the vectorization label 
subscripts unless referring to a general cr-basis. The main 
transformation we will be interested in is then from col- 
vec to another arbitrary orthononormal operator basis 
{do,}. Tensor networks for the change of basis T c ^ a and 
its inverse T a ^, c are respectively illustrated in Fig. |21| 




(a) Col-vec to cr-basis 




(b) cr-basis to col-vec 

FIG. 21. Vectorization change of basis operators between col- 
vec convention and an orthonormal basis {cr a } for £.(X,y). 

As was illustrated in Fig. [l7j in the case where the 
second basis is the row-vec convention we have 



= SWAP. 



(20) 



One final important result that often arises when deal- 
ing with vectorized matrices is Roth's Lemma for the 
vectorization of the matrix product ABC [T3]. Given 
matrices A,B,C£ C(X) we have 



| ABC)) = {C T ®A)\B)) 



(21) 



The graphical tensor network proof of this lemma is 
found in Fig. [22} 



III. REPRESENTATIONS OF COMPLETELY 
POSITIVE MAPS 

Before presenting the main results found in this study 
we first recall several common mathematical descriptions 
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ABC 



FIG. 22. Graphical proof of Roth's lemma (21 1 in col-vec 
convention. 



for completely-positive trace-preserving maps, and show 
how several key properties may be captured graphically 
using the diagrammatic notation we have developed spe- 
cific for this purpose. The representations we will con- 
sider are the Kraus (or operator-sum) representation, the 
system-environment (or Stinespring) representation, the 
Liouville superoperator description based on vectoriza- 
tion of matrices, and the Choi-matrix or dynamical ma- 
trix description based on the Choi-Jamiolkowski isomor- 
phism. We will also describe the often used process ma- 
trix representation and show how this can be considered 
as a change of basis of the Choi-matrix. Following this, in 
Section Hvl we will show how our framework now enables 
one to freely transform between these representations as 
illustrated in Fig. [T] 



A. Kraus / Operator-Sum Representation 

The first representation of CPTP-maps we cast into 
our framework is the Kraus |16j or operator-sum pQ rep- 
resentation. This representation is particularly useful in 
phenomenological models of noise in quantum systems. 
Kraus's theorem states that a linear map £ £ T(X, y) is 
CPTP if and only if it may be written in the form 



(22) 



where the Kraus operators {K a : a = 0, D— 1}, K a € 
C(X,y), satisfy the completeness relation 



D 



(23) 



a=0 



This Kraus evolution has the graphical representation 
shown in Fig. [23] 



6(P) 



= — K - p - K< — 



FIG. 23. Tensor network for the evolution of a state p under 
a CP-map £ in the Kraus representation (|22|). 



The maximum number of Kraus operators needed for 
a Kraus description of £ is equal to the dimension of 
C{X,y). For the case where X = y = C d the maximum 



number of Kraus operators is d 2 , and the minimum num- 
ber case corresponds to unitary evolution where there is 
only a single Kraus operator. 

It is important to note that the Kraus representation of 
£ is not unique as there is unitary freedom in choosing the 
Kraus operators. We can give preference to a particular 
representation called the Canonical Kraus Representa- 
tion |17j which is the unique set of Kraus operators satis- 
fying the orthogonality relation Tr[K^K^] = X a S a /3. The 
canonical Kraus representation will be important when 
transforming between representations in Section |IV| 



B. System-Environment / Stinespring 
Representation 

The second representation of CPTP-maps we consider 
is the system-environment model pQ, which is typically 
considered the most physically intuitive description of 
open system evolution. This representation is closely re- 
lated to (and sometimes referred to as) the Stinespring 
representation as it can be thought of as an application 
of the Stinespring dilation theorem 18J, which we also 
describe in this section. In this model, we consider a sys- 
tem of interest X, called the principle system, coupled 
to an additional system Z called the environment. The 
composite system of the principle system and environ- 
ment is then treated as a closed quantum system which 
evolves unitarily. We recover the reduced dynamics on 
the principle system by performing a partial trace over 
the environment. Suppose the initial state of our compos- 
ite system is given by /)®t e C{X®Z), where r G C{Z) 
is the initial state of the environment. The joint evolu- 
tion is described by a unitary operator U & C(X ® Z) 
and the reduced evolution of the principle system's state 
p is given by 



£{p) = Ti z [U{p®T)tf] 



(24) 



For convenience we can assume that the environment 
starts in a pure state r = |wo)(i | o|i and in practice one 
only need consider the case where the Hilbert space de- 
scribing the environment has at most dimension d 2 for 
X = C d pQ. The system-environment representation of 
the CP-map £ may then be represented graphically as 



shown in Fig. 24 



- Sip) - = 




FIG. 24. Tensor network for the evolution of a state p by 
a CP-map £ in the system-environment representation ( |24[ |. 
The environment is initially in the state \vq) and the joint 
evolution is given by a unitary operator U. 
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The system-environment model is advantageous when 
modelling the environment as a physical system. How- 
ever, care must be taken when ascribing physical reality 
to any particular model as the system-environment de- 
scription is not unique. This is not surprising as many 
different physical interactions could give rise to the same 
reduced dynamics on the principle system. This free- 
dom manifests in an ability to choose the initial state 
of the environment in the representation and then ad- 
just the unitary operator accordingly. In practice, the 
system-environment model can be cumbersome for per- 
forming many calculations where the explicit dynamics 
of the environment system are irrelevant. The remaining 
descriptions, which we cast into diagrammatic form, may 
be more convenient in these contexts. 

Note that the system-environment evolution for the 
most general case will be an isometry and this is cap- 
tured in Stinespring's representation 18J. Stinespring's 
dilation theorem states that a CP-map £ € C(X,y) can 
be written in the form 



£(p) = Tr z [Aptf] 



(25) 



where A G C(X,y<E)Z) and the Hilbert space Z has 
dimension at most equal to C(X,y). Further, the map 
£ is trace preserving if and only if A^A = \x |18j . 

In the case where y = X, the Hilbert space X ® Z 
mapped into by the Stinespring operator A is equivalent 
to the joint system-environment space in the system- 
environment representation. Hence, as illustrated in 
Fig. |25| one may move from the system-environment de- 
scription to the Stinespring representation by letting 

A = U -{\ x ®\v Q )\ (26) 

where Voq) £ Z is the initial state of the environment. 



-S(p)-= i 







p 


u 







Iff 



A 9 A 1 



FIG. 25. Graphical depiction of the relationship between the 



Stinespring ( 25 \ and system-environment ( 24 1 representations 



of a CP map £. 

This close relationship is why these two representa- 
tions are often referred to by the same name, and as we 
will show in Section |IV E[ it is however straight forward 
to construct a Stinespring representation from the Kraus 
representation. However, generating a full description of 



the joint system-environment unitary operator U form 
a Stinespring operator A is cumbersome. It involves an 
algorithmic completion of the matrix elements in the uni- 
tary U not contained within the subspace of the initial 
state of the environment [T7] . Since it usually suffices to 
define the action of U when restricted to the initial state 



of the environment, which by Eqn (26) is the Stinepsring 
representation, this is often the only transformation one 
need consider. 

A further important point is that the evolution of the 
principle system £ (p) is guaranteed to be CP if and only 
if the initial state of the system and environment is sep- 
arable; pxz — Px ® Pz- In the case where the physi- 
cal system is initially correlated with the environment, 
it is possible to have reduced dynamics which are non- 
completely positive |19H21| . however such situations are 
beyond the scope of this paper. 



C. Louiville-Superoperator Representation 

We now move to the linear superoperator or Liouville 
representation of a CP-map £ € C(X,y). This represen- 
tation is based on the vectorization of the density matrix 
p h> \p)) with respect to some orthonormal ope rator 
.,d 2 



II B 



basis {a a : a = 0, ...,d 2 ~ 1} introduced in Section 
Once we have chosen a vectorization basis (eol-vee in our 
case) we define the superoperator for a map £ £ T(X, y) 
to be the linear map 

S : X <gi X -> y <g) y : \p)) h> \£{p))) (27) 

This is illustrated in Fig. [26} 




FIG. 26. Tensor network for the evolution of a s tate p by 
CP-map £ in the superoperator representation ( 27 1 . 



In the eol-vee basis we can express the evolution of a 
state p in terms of tensor components of S as 



(28) 



For the case where £ € T(X), it is sometimes use- 
ful to change the basis of our superoperators from the 
eol-vee basis to an orthonormal operator basis {<J a } for 
C{X). This is done using the basis transformation oper- 
ator Tc-j-cr introduced in Section |il B| We have 



S ■ Tj_ 



J2 Sa P l Cr ">)(( Cr /3| 

a/3 



(29) 
(30) 



where the subscript a indicates that S a is the superop- 
erator in the cr-vec convention. The tensor networks for 
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this transformation is given in Fig. 27 Note that for 
a general map £ G T(X,y) we could do a similar con- 
struction but would need different bases for the initial 
and final Hilbert spaces C(X) and C(y). 





















So- 




T 




s 




ft 



















FIG. 27. Tensor network for changing vectorization conven- 
tion of a superoperator from col-vec to an arbitrary operator 
basis {a a }. 

The structural properties the superoperator S must 
satisfy for the linear map £ to be trace-preserving (TP), 
hermitian-preserving (HP) and completely positive (CP) 
are [T7j: 



£ is HP 



S = S k 



(31) 



£ is TP S. 



_ s p = xi s IX (32) 

(33) 



mm.nv u n7y 



£ is CP 



(34) 
(35) 



> <Szi8>£|/=>Ab)) > VpAB > 

Note that there is not a convenient structural criteria 
on the superoperator S which specifies if £ is a CP-map 
to test for positivity or complete positivity one generally 
uses the closely related Choi Matrix representation. 

Superoperators are convenient to use for many prac- 
tical calculations. Unlike the system-environment model 
the superoperator S is unique with respect to our choice 
of vectorization basis. Choosing an appropriate basis to 
express the superoperator in can often expose certain in- 
formation about a quantum system. For example, if we 
want to model correlated noise for a mutli-partite system 
we can vectorize with respect to the mutli-qubit Pauli 
basis. Correlated noise would then manifest as non-zero 
entries in the superoperator corresponding to terms such 
as o~ x ®o~ x . One can perform this transformation by using 
the tensor network shown in Fig. |27| 



For X = C d , the explicit construction of the Choi- 
matrix is given by 



d-l 



Ar 



£ \i){j\®£(\i){j\) 



A r = 



d-l 

E 

i,j=0 



£(\i)(j\)®\*)(3\ 



(38) 



(39) 



where {\i} : i = 0, . . . , d — 1} is an orthonormal basis for 
X. 

We call the two conventions col-A and row-A due to 
their relationship with the vectorization conventions in- 
troduced in Section III Bl The Choi-Jamiolkowski iso- 
morphism can also be thought of as having a map 
£ G T(X,y) act on one half of a Bell state |$+) = 
J2i K) ® I*) e ® X), and hence these conventions 
corresponding to which half of the Bell state it acts on: 

A c = (X®£)\$+)($+\ (40) 

A r = (£ ®Z)|$ + ><$+| (41) 

where X € T(X , X) is the identity map. In what follows 
we will use the col-A convention and drop the subscript 
from A c . We note that the alternative row-A Choi-matrix 
is naturally obtained by applying the bipartite-SWAP 
operation to A c . 

As will be considered in Section IV C[ if the evolution 
of the CP map £ is described by a Kraus representation 
{Ki}, then the Choi-Jamiolkowski isomorphism states 
that we construct the Choi-matrix by acting on one half 
of a bell state with the Kraus map. This is illustrated in 
Fig. [281 





D. Choi-Matrix Representation 

The final representation shown in Fig. [l] is the Choi 
matrix |22j . or dynamical matrix [17] . This is an ap- 
plication of the Choi-Jamiolkowski isomorphism which 
gives a bijection between linear maps and linear opera- 
tors [23j . Similarly to how vectorization mapped linear 
operators in £(X, y) to vectors in X ® y or y ® X, the 
Choi-Jamiolkowski isomorphism maps linear operators in 
T(X,y) to linear operators in C(X <g> y) or £(y ® X). 
The two conventions are 



col-A :T(X,y) 
row-A : T{X,y) 



C(X(g>y) : £ ^ A c (36) 
C(y®X): £^ A r . (37) 



FIG. 28. A graphical depiction of the Choi-Jamiolkowski iso- 
morphism for a Kraus description of a CP map £ . Note that 
in general any tensor network describing a linear map 8, not 
just the Kraus description, may be contracted with one-half 
of the maximally entangled state |$ + }($ + | to construct the 
Choi-matrix. 

With the Choi-Jamiolkowski isomorphism defined, the 
evolution of a quantum state in terms of the Choi-matrix 
is then given by 

£{p) = Tr x [(p T ® l y )A] (42) 

or in terms of tensor components 



£(j>)t 



= > A 



n.m 



(43) 
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The tensor network for Eqn (42 1 is given in Fig. [29 



and the graphical proof for the case where £ is described 
by a Kraus representation is given in Fig. [30} 



-~B(p)\- = Cnr ttl) = CrfV 



A 



A 



FIG. 29. Tensor network for the evolution of a state p by 
CP-map £ in the Choi-matrix representation (42 I. 



A 



K 



Due to the similarity of vectorization and the Choi- 
Jamiolkowski isomorphism, one could then ask what hap- 
pens if we vectorize in a different basis. This change of 
basis of the Choi-matrix is more commonly known as the 
X-matrix which we will discuss next. However, such a 
change of basis does not change the eigen-spectrum of a 
matrix, so the positivity criteria in Eqn ( 48 1 holds for 
any basis. 

Another desirable property of Choi matrices is that 
they can be directly determined for a given system 
experimentally by ancilla assisted process tomography 
(AAPT) [2~4l, I25j . This is an experimental realization of 
the Choi-Jamiolkowski isomorphism where one prepares 
the maximally entangled state A J2i=o I*) ® N) on tnc 
joint ancilla-system and subjects one half to the channel 
£. State tomography is then performed on the output 
and the reconstructed state is equal to the Choi matrix 
up to the normalization factor d. 



— H" p h 



6(p) 



FIG. 30. Proof of 
CP-map £ (Fig. [23 1 



)f ( 42 1 assuming a Kraus representation of a 
[23 1 and using the snake-equations (Fig. 



The structural properties the Choi-matrix A must sat- 
isfy for the linear map £ to be trace-preserving (TP), 
hermitian-preserving (HP) and completely positive (CP) 
are \TT\: 



£ is HP 



A T = A 



A 



A 



£ is TP Try [A] = 1. 



A 



£ is CP <^=> £ is CP 



A > 0. 



(44) 



(45) 



(46) 



(47) 



(48) 



The Choi-matrix for a given map £ is unique with re- 
spect to the isomorphism convention chosen. We will pro- 
vide tensor networks to illustrate a close relationship to 
the superoperator formed with the corresponding vector- 
ization convention in Section FlV A[ The Choi-matrix finds 
practical utility as one can check the complete-positivity 
of the map £ by computing the eigenvalues of A, and it 
is also necessary to construct the Choi-matrix for a given 
superoperator to transform to the other representations. 



E. Process Matrix Representation 

As previously mentioned, one could consider a change 
of basis of the Choi-matrix analogous to that for the su- 
peroperator. The resulting operator is more commonly 
known as the \-matrix or process matrix pQ. Consider 
Hilbert spaces X = C dx , y = C d « and let D = d x d y . 
If one chooses an orthonormal operator basis {a a : a = 
0, D — 1} for C(X, y), then a CPTP map £ G C(X, y) 
may be expressed in terms of a matrix % G ® y) as 



D-l 



(49) 



a,P=0 



where the process matrix \ is unique with respect to the 
choice of basis {er Q }. 

In terms of the Choi-matrix, we have that the process 
matrix with respect to an orthonormal operator basis 
{a a } is given by 



= T* -A -Tr- 



et, 



a 



(50) 
(51) 



where T c ^ a is the vectorization change of basis operator 
introduced in Section [ITB] Thus evolution in terms of the 
X-matrix is analogous to our Choi evolution as given in 

Fig. Em 



The graphical proof asserting the validity of Eqn ( 50 1 is 



given in Fig. [32j We also see that if one forms the process 
matrix with respect to the eol-vee basis a a = Ejj where 
a = i + dj and d is the dimension of Ji, then x c = A. 

Since the process matrix is a unitary transformation of 
the Choi-matrix, it shares the same structural conditions 
for hcrmitian preservation and complete-positivity as for 
the Choi-matrix given in Eqns ( 44 ) and ( 48 ) respectively. 



The condition for it to be trace preserving is different 
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— £(p) — = 





- P - 




* T 


X 


- A c-«r . 



where the superscripts a, uj denote the basis of the x~ 



matries. This is illustrated in Fig. 33 

















x w 






x a 




T 



FIG. 31. Tensor network for the evolution of a state p by CP- 
map £ in the ^-matrix representation 1 49 1 . It is analogous to 



evolution in the Choi-matrix representation (Fig. 29 1 



FIG. 33. Tensor network for changing a process matrix x' 7 
with respect to basis {u a } to a process matrix x" with respect 
to basis {uj a }. 



— 6(p) — = — cr - p - a* — 



X 









X 








X 





A 



FIG. 32. Graphical proof of the equivalence of the Choi- 
matrix and process matrix under change of basis given in 
(501. 



however. These conditions are 



£ is TP 
£ is HP 
£ is CP 



Try [T e ^ a x T Ua] 
X 1 = X 

x>o. 



(52) 

(53) 
(54) 



To convert a process-matrix % in a basis {cr a } to an- 
other orthonormal operator basis {w Q }, we may use the 
inverse transformation to the superoperator change of ba- 
sis from Section llH CI That is 



T 1 ' • v CT • T _^ 

a/3 



(55) 
(56) 



IV. TRANSFORMING BETWEEN 
REPRESENTATIONS 

We now proceed to the task of describing how one may 
transform between representations of completely-positive 
trace-preserving maps, as depicted in Fig. [T] In particu- 
lar, the transformations depicted as solid arrows in Fig.[T] 
have succinct descriptions in the graphical calculus we in- 
troduced in Section [III These transformations are based 
on the wire bending dualities for reshuffling, vectoriza- 
tion, and the Choi-Jamiolkowski isomorphism. While 
the remaining transformations depicted as dashed lines 
are not based on dualities, but rather non-linear decom- 
positions, or constructions, they also have diagrammatic 
representations in our graphical calculus for completely 
positive maps. 



A. Transformations between the Choi-matrix and 
superoperator representations 

The Choi-matrix and superoperator are naturally 
equivalent under the reshuffling wire bending duality in- 
troduced in Section II A In the col (row) convention we 
may transform between the two by applying the bipar- 
tite col (row)-reshuffling operation R introduced in Sec- 
Let A € C(X (g) y) be the Choi-matrix, and 
® Xy (g) y) be the superoperator, for a map 



tion 



II A 



S eC[X 
£ e T{X,y). Then we have 

A = S R 

This is illustrated in Fig. 34 
nents we have 



S = A 1 



(57) 

In terms of tensor compo- 



(58) 



where m, n and ^, v index the standard bases for X and 
y respectively. 

Graphical proofs of the relations A Rc — S and S Rc = A 



are given in Fig.s 35 and 36 respectively. 

To transfer between a process matrix with respect to 
an arbitrary operator basis, and a superoperator with 
respect to an arbitrary vectorization basis, we must first 
convert both to col-vec (or row-vec) convention and then 
proceed by reshuffling. 
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A 



S 



S 



A 



FIG. 34. Transforations between superoperator and Choi ma- 
trix are done by the wire bending duality for the bipartite 
reshuffling operation (Fig. 15f|. 



A 



A 
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G 
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G 
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6(P) 



G^O 



A 



FIG. 36. Graphical proof that the reshuffled superoperator 
S R satisfies the evolution equation for the Choi-matrix rep- 
resentation (|42|). 




FIG. 35. Graphical proof that the reshuffled Choi-matrix A R 
satisfies the evolution equation for the superoperator repre- 
sentation (271. 



Note that reshuffling is it's own inverse, ie (A R ) R = A, 
hence the solid bi-directional arrow connecting the Choi- 
matrix and superoperator representations in Fig. [l] This 
is the only transformation between the representations 
we consider which is linear, bijective, and self-inverse. 



B. Transformations to the superoperator 
representation 

Transformations to the superoperator from the Kraus 
and system-environment representations of a CP-map are 
also accomplished by a wire-bending duality, in this case 
vectorization. However, unlike the bijective equivalence 
of the Choi-matrix and superoperator under the reshuf- 
fling duality, the vectorization duality is only surjective. 

If we start with a Kraus representation for a CPTP 
map E £ C(X,y) given by {K a : a = 0, ...,£> - 1}, 
with K a € C(X, y), we can construct the superoperator 

s eC{x®x,y®y)hy 



D-l 



(59) 



The corresponding tensor network is given in Fig. [37] and 
the graphical proof of this relationship follows directly 
from Roth's lemma and is found in Fig. 38 
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K 



K 



FIG. 37. Tensor network for constructing the superoperator 
representation S for a CPTP-map £ from a Kraus represen- 
tation {K a : a = 0, . . . D — 1}. 




K - p - K 




FIG. 38. Graphical proof asserting the validity of Eqn ( 59 \ for 



constructing a superoperator representation S from a Kraus 
representation {K a }. 



Starting with a system-environment (or Stinespring) 
representation with a map £ G C(X i y) with system 
Hilbert space X = C d and environment Hilbert space 
Z = C D with 1 < D < dim(C(X ,y)), we may construct 
the superoperator for this map from the joint system- 
environment unitary U and initial environment state \vq) 
by 



S = Y,{a\U\v )® {a\U\v ) 



(60) 



where {\a) : a = 0, D — 1} is an orthonormal basis for 
Z. The corresponding tensor network is given in Fig. |39| 
As with the Kraus to superoperator transformation, the 



proof of Eqn ( 60 1 follows from Roth's lemma. 




U 



U 



FIG. 39. Tensor network for constructing the superoper- 
ator representation 5 for a CPTP-map £ from a System- 
Environment representation with initial environment state 
| no) and joint unitary U. 



the superoperator is unique, while both the Kraus and 
system-environment representations are not. Hence we 
have solid single directional arrows in Fig. [T] connecting 
both the Kraus and system-environment representations 
to the superoperator. The inverse transformation from a 
superoperator to the Kraus or system-environment rep- 
resentation requires a canonical decomposition of the op- 
erator S (via first reshuffling to the Choi-matrix), which 
is detailed in Sections IIVDI and IIV El 



C. Transformations to the Choi-matrix 
representation 

Transforming to the Choi-matrix from the Kraus 
and system-environment representations is accomplished 
via a wire-bending duality which captures the Choi- 
Jamiolkowski isomorphism. As with the case of trans- 
forming to the superoperator, this duality transformation 
is surjective but not injective. 

Given a set of Kraus matrices {K a : a = 0, D — 1} 
where K a e C{X ,y) for a CPTP-map £ € C(X,y), one 
may form the Choi-Matrix A as was previously illustrated 
in Fig.[28]in Section[lII D| In terms of both Dirac notation 



and tensor components we have: 



a=e i\i){j\®j2 K ^M K i) ( 61 ) 

i.j \ a / 

= Y / \K a ))((K a \ (62) 

a 

^(^ Q ) Mm (^ Q ) ra . (63) 



where {\i}} is an orthonormal basis for X, m, n index the 
standard basis for X, and fi, v index the standard basis 

for y. 

Given a system-environment representation with joint 
unitary U £ C(X , Z) and initial environment state |^o) € 
Z we have 

A = (1*>0'I ® Tr ^ [Um\ ® MMU^}) (64) 
This is shown graphically in Fig. |40] 




U 







Note that while the vectorization wire bending duality 
is invertible, these transformations to the superoperator 
from the Kraus and system-environment representations 
are single directional. In both cases injectivity fails as 



FIG. 40. Construction of the Choi-matrix representation 
A from a system-environment representation with initial 
environment state \v) and joint unitary U via the Choi- 
Jamiolkowski isomorphism 1401). 
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The proof of these transformations follow directly 
from the definition of the Choi- matrix in Eqn ( |38[ ), and 
the tensor networks for the evolution via the Kraus or 
system-environment representations given in Figues [23] 



and 24 respectively. As with the vectorization transfor- 
mation to the superoperator discussed in Section |IVB[ 
even though the Choi-Jamiolkowski isomorphism is lin- 
ear these transformations are single directional as injec- 
tivity fails due to the non-uniqueness of both the Kraus 
and system-environment representations. Hence we have 
the solid single-directional arrows connecting both the 
Kraus and system-environment representations to the 
Choi-matrix in Fig. [I] 

This completes our description of the linear trans- 
formations between the representations of CP-maps in 
Fig. [TJ We will now detail the non-linear transformations 
to the Kraus and system environment representations. 



D. Transformations to the Kraus Representation 

We may construct a Kraus representations from the 
Choi-matrix or system-environment representation by 
the non-linear operations of spectral-decomposition and 
partial trace decomposition respectively. To construct a 
Kraus representation from the Superoperator however, 
we must first reshuffle to the Choi-matrix. 

To construct Kraus matrices from a Choi matrix we 
first recall the graphical Spectral decomposition we intro- 
duced as an example of our colour summation convention 
in Section [II] If £ is CP by Eqn Q we have A £ > and 
hence the spectral decomposition of the Choi-matrix is 
given by 




L J 



A 




P N At 



L _ _ J L _l 



K - p - 10 



FIG. 42. Graphical proof of Fig. [41] for the construction of 
a Kraus representation {K a } for a CPTP-map £ from the 
Choi-matrix representation A. 



(65) 



where X a > 0. Hence we can define Kraus operators 
Ka = A Q A Q where A a is the unique operator satisfy- 
ing \A a )) = \(f> a ) as illustrated in Fig. 
of Kraus operators will be equal to t 



41 The number 
le rank r of the 



Choi-matrix, where 1 < r < dim(£(A', y)). The graphi- 
cal proof of this transformation is given in Fig. [42] The 
proof that Kraus operators satisfy the completeness re- 
lation follows from the trace preserving property of A in 
Eqn ((io). This is given in Fig. [43 



G 




Note that since A, and the ^-matrix are related by a 
unitary change of basis, the Kraus representations con- 
structed from their respective spectral decompositions 
will also be related by the same transformation. Each will 
give a unitarily equivalent Canonical Kraus representa- 
tion of £ since the eigen-vectors are orthogonal. Thus we 
have described the arrow in Fig. [I] connecting the Choi- 
matrix to the Kraus representation. It is represented as 
a dashed arrow as it involves a non-linear decomposition, 
and is single directional as this representation transfor- 
mation is injective, but not surjective. Surjectivity fails 
as we can only construct the canonical Kraus represen- 
tations for £ . The reverse transformation is given by the 
Jamiolkowski isomorphism described in Section ]"lV C| 



FIG. 41. Construction of a Kraus representation K a for a 
CPTP-map £ from the Choi-matrix representation A by the 
de-vectorization of the spectral decomposition of the Choi- 
matrix A = Y, a ^l\<Pa){4> a \- 



Starting with a system-environment representation 
with joint unitary U G C{X , Z) and initial environment 
state \vq) S -Z, we first choose an orthonormal basis 
{\a) : a = 0,...,D — 1} for Z. We then construct the 
Kraus representation by decomposing the partial trace 
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— f? - K — = 





<t>X * 




K a 




U 






4 







FIG. 44. Construction of a Kraus representation {Kj : a = 
0, ...,£> — 1} from a system-environment representation with 
initial environment state \vo), joint unitary U, and an or- 
thonormal basis for the environment {|a} : a = 0, ... ,D — 1}. 




c 



J 



FIG. 43. Graphical proof of the completeness relation ( 23 1 for 
the Kraus representation {K a } constructed from the Choi- 
matrix representation as given in Fig. |41| This follows from 
the trace preserving property of the Choi-matrix representa- 
tion pfk. 



in this basis as follows 



£(p) = Tr E [U(j>®\v)(v\)rf] 

D-l 

= {a\U\v )p(v \U^\a) 



a=0 
D-l 



= K aP Kl 



(67) 



(68) 



A sequence of tensor contractions illustrating the validity 
of this results are given in Fig. |45| Hence we may define 
Kraus matrices 
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FIG. 45. Graphical proof Proof of Fig. [44] for constructing a 
Kraus representation {K a } from a system environment repre- 
sentation with initial environment state |t>o) and joint unitary 
U. 



(a\U\vo) 



(69) 



leading to the tensor network given in Fig. |44j 

Though the Kraus and system-environment represen- 
tations are both non-unique, for a fixed basis for the en- 
vironment this partial trace decomposition is an injec- 
tive transformation between the Kraus and Stinespring 
representations (or equivalently between the Kraus and 
system-environment representations when the joint uni- 
tary is restricted to a fixed initial state of the environ- 
ment). To see let {K a } and {J a } be two Kraus repre- 
sentations constructed from Stinespring representations 



A and B respectively for a CPTP-map £ e C(X,y). 
Then we have that 



K — 1 

^* Ai a j B} a j 

&A = B. 



(70) 
(71) 
(72) 
(73) 



Since the Stinespring operators satisfy A = U\vq) and 
B = V\vq) for some joint unitaries U and V, we must 
have that Uq — Vq where Uq and Vq are the joint unitaries 
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restricted to the subspace of the environment spanned by 
|«o>. 

This transformation can be thought of as the re- 
verse application of the Stinespring dilation theorem, and 
hence for a fixed choice of basis (and initial state of the 
environment) it is invertible. The inverse transforma- 
tion is the Stinespring dilation, and as we will show in 
Section [IV E| since the inverse transformation is also in- 
jective this transformation is a bijection. However, since 
the partial trace decomposition involves a choice of basis 
for the environment it is non-linear — hence we use a 
dashed bi-directional arrow to represent the transforma- 
tion from the system-environment representation to the 
Kraus representation in Fig. [I] 



E. Transformations to the system-environment 
representation 

We now describe the final remaining transforma- 
tion given in Fig. [I] the bijective non-linear transfor- 
mation from the Kraus representation to the system- 
environment, or Stinespring, representation. The 
system-environment representation is the most cumber- 
some to transform to since it is it involves the unitary 
completion of a Stinespring dilation of a Kraus represen- 
tation. Thus from a Superoperator one must first reshuf- 
fle to the Choi-matrix, and form the Choi-matrix descrip- 
tion one must then spectral decompose to the canoni- 
cal Kraus representation before finally constructing the 
system-environment as follows. 

Let {K a : a = 0, ...,D - 1} where 1 < D < 
dim(C(X,y)), be a Kraus representation for the CP- 
map £ £ T(X,y). Consider an ancilla Hilbert space 
Z = C D , this will model the environment. If we choose 
an orthonormal basis for the environment, {\a) : a = 
0, ...,£> — 1}, then by Stinesprings dilation theorem we 
may construct the Stinespring matrix for the CP map £ 
by 



D-l 

a=Q 



® \a) 



(74) 



Recall from Section |III B| that the Stinespring repre- 
sentation is essentially the system-environment represen- 
tation when the joint unitary operator is restricted to the 
subspace spanned by the initial state of the environment. 
Hence if we let \vq) £ He be the initial state of the en- 
vironment system, then this restricted unitary is given 
by 



U a = J2 Ka ® \ a )( v o\ 



cumbersome and is unnecessary to describe the evolution 
of the CP-map £ [17]. 





K 



(a) Stinespring operator 





Ya 



(b) Restricted unitary 



FIG. 46. Construction of a Stinespring representation A and 
system-environment representation with initial environment 
state |i>o) and joint unitary U for a CPTP-map £ from a 
Kraus representation {Ka}. Note that Uo is the restriction 
of the joint unitary U to the subspace of the environment 
spanned by the initial state \vo). 

Hence we have finished characterizing the final trans- 
formations depicted in Fig. [I] connecting the Kraus rep- 
resentation to the system-environment representation by 
Stinespring dilation. As previously mentioned in Sec- 
tion IV D[ for a fixed choice of basis for the environment, 
and initial state of the environment, the transformation 
between Kraus and Stinespring representations is bijec- 
tive (and hence so is the transformation between Kraus 
and system-environment representations when restricted 
to the subspace spanned by the initial state of the envi- 
ronment). Though both these representations are non- 
unique, by fixing a basis and initial state for the environ- 
ment we ensure that this transformation is injective. To 
see this let Uq and Vo be unitaries restricted to the state 
\vo) constructed from Kraus representations, {K a } and 
{J a } respectively, for £ £ C(X,y). Then 

a a 



<&Kp = J, 







(78) 



(75) 



Bijectivity then follows from the injectivity of the inverse 
transformation — the previously given construction of a 
Kraus representation by the partial trace decomposition 
of a joint unitary operator in Fig. |44| 



as illustrated in Fig. [46] The graphical proof of this con- 
struction giving the required evolution of a state p can be 
found in Fig. |47[ In principle, one may complete the re- 
maining entries of this matrix to construct the full matrix 
description for the unitary U, however such a process is 



V. CONCLUSION 

The study of completely-positive trace-preserving 
maps is an old topic, so it is perhaps surprising that 



18 











■ p - 




u. 


• 
















K - p - 



IP 



K - f> - fC 





Sip) 









FIG. 47. Graphical proof that the restricted unitary Uo con- 
structed from a Kraus representation by ( 75 1 gives the correct 
description of CP map £ . 



there are still new insights to be gained by investigat- 
ing their structure using new techniques. Further, while 
the application of CPTP-maps to describing the evolu- 
tion of open quantum systems is well understood, it is 
a surprisingly difficult task to find a concise summary 
of the properties of, and transformations between, their 
various mathematically equivalent representations used 
in the quantum information processing literature. The 
graphical calculus for open quantum systems presented 
in this paper has enabled us to unify, and hence trans- 
form freely between, the various common representations 
of CPTP-maps by performing diagrammatic manipula- 
tions of their respective tensor networks. A summary of 
these transformations between the different representa- 
tions was given in Fig. [T] 

We found that many of these transformations between 
representations of CPTP-maps corresponded to wire 
bending dualities in our graphical calculus, which have a 
particularly succinct tensor network description. These 
transformations are depicted by solid arrows between two 
boxes labelling representations in Fig. [T] Of these dual- 
ity transformations, only the reshuffling operation con- 



necting the Choi-matrix and Liouville-superoperator is 
bi-directional — the reshuffling operation is bijective and 
self-inverse, and hence the same transformation takes the 
Choi-Matrix to the superoperator as takes the superop- 
erator to the Choi-matrix. The two other wire bend- 
ing dualities are vectorization, which transforms both the 
Kraus and system-environment representations to the su- 
peroperator representation, and the Choi-Jamiolkowski 
isomorphism, which transforms the same two representa- 
tions to the Choi-matrix. These duality transformations 
are only single directional as they are not injective, and 
hence depicted by a one-way arrow connecting the ap- 
propriate boxes labelling these representations in Fig. [T] 
The reason these transformations are single directional, 
as apposed to the bi-directional transformation between 
the Choi-matrix and superoperator, is due to the non- 
uniqueness of the Kraus and system-environment rep- 
resentations of a CPTP-map. The transformation is a 
many-to-one (surjective) mapping and not strictly invert- 
ible without first specifying some form of decomposition 
of the superoperator or Choi-matrix. 

The transformations we presented for transitioning 
from the Choi-matrix to the Kraus representation, and 
between the Kraus and system-environment representa- 
tions, were not based solely on wire bending dualities. 
These transformations are depicted by dashed arrows in 
Fig. [T] where the dash is meant to indicate that they 
are non-linear transformations. This non-linearity arose 
from the decompositions and constructions involved, for 
example the spectral decomposition of a positive-semi 
definite operator in the Choi-matrix to Kraus representa- 
tion transformation. In our case, these non-linear trans- 
formations were all also one directional due to the non- 
uniqueness of the representation being transformed to. 
There is unitary freedom in constructing them — for the 
Choi- matrix to Kraus representation transformation, one 
could change the basis of the eigenvectors with respect 
to a vectorization convention and still arrive at a valid 
Kraus representation; for Kraus to the system environ- 
ment representation one may choose any orthonormal ba- 
sis in the construction of the joint system-environment 
unitary in Eqn ( 75 1 ; and for the system environment 



to Kraus representation one may decompose the partial 
trace over the environment in any orthonormal basis. 

Having new tools to investigate old problems can of- 
ten lead to surprising new results, and we believe there 
are many potential applications in QIP for the graphical 
calculus we have presented in this paper. 
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Appendix A: Tensor networks 

We will now prove the consistency of several of the 
basic tensor networks introduced in Section [XTJ and in 
doing so illustrate how one may use our graphical calculus 
for diagrammatic reasoning. 

The colour summation convention we have presented 
represents diagrammatic summation over a tensor index 



by colouring the appropriate tensors in the diagram. In 
this convention summation over a Kronecker delta is as 
shown in Fig. |48| This expression is used in several of 
the following proofs. 



O <H> 



FIG. 48. Colour summation convention for ^2% j = 
. Sij. Multiplication of a tensor network by this quantity 
converts all diagrams with either of the two colours present 
in the left hand side, to the same colour. 



We begin with the proof of the trace of an operator A 



which is given in Fig. 49 



VT7V = <JTl-> 



= 4uU> 



FIG. 49. Proof that our tensor network for graphical trace 
does return Tr(A) = ^ An. 

For illustrative purposes, to prove this algebraically we 
note that the tensor networks for trace correspond to the 
algebraic expressions ($ + |^4® and ($ + |l<g)^4|$ + ), 

and that 



(<t>+\A®l\$+) = J2(i\A\j)(i\j) 

i,j 

i 

= Tr[A\. 



(Al) 

(A2) 

(A3) 
(A4) 



Similarly we get ($+|l ® A\$+) = Tr[A]. 

To prove the snake equation we must first make the 
following equivalence for tensor products of the elements 
\i) and 



01®l*) = l*)®0'l = l*>01 



(A5) 



This is illustrated diagrammatically in Fig. [50] 

With this equivalence made, the proof of the snake- 
equation for the "S" bend is given in Fig. [5l] The 
proof for the reflected "S" snake-equation follows nat- 
urally from the equivalence defined in Fig. |50| 
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FIG. 50. Equivalence of tensor product expressions \i) 
01, 01 ® K>i l*>01- 
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FIG. 51. Graphical proof of the snake equation. 
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FIG. 52. Proof of transposition tensor network. 



The proof of our tensor network for the transposition 
of a linear operator A is shown in Fig. |52| To prove this 
algebraically we note that the corresponding algebraic 
equation for the transposition tensor network is 



= 1® ($+|(l(g) A ® <i 


5 1 (A6) 


id 


(A7) 


= Y^{j\Am){j\ 

id 


(A8) 


= J2{i\A T \j)\i){j\ 

id 


(A9) 


= j2m\A T \j)(j\ 

id 


(A10) 


= A T . 


(All) 



The proof for transposition by counter-clockwise 
wire bending follows from the equivalence relation in 
Eqn dAlb. 



With the tensor network for transposition of an op- 
erator proven, the proof of transposition by contracting 
through a Bell-state |$ + ) is then an application of the 



snake-equation as shown in Fig. 53 



— A T 



A 







FIG. 53. Proof of transposition by "sliding" . 

Appendix B: Vectorization change of basis 

We now prove that the vectorization change of basis 
operator T a ^ u indeed functions as claimed. Let H be a 
e?-dimensional complex Hilbert space and let {er Q : a = 
0,...,d 2 — \},{uj a : a = 0,...,d 2 — 1} be orthonormal 
operator bases for C(H). Define an operator T a ^ u £ 



(Bl) 



where {\a) : a = 0, d 2 — 1} is the computational basis for 
H®H. 
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We claim that for any linear operator A e The proof is as follows: 



T a ^\A)) a =(^2\*)((^ a \)j\A)) a (B3) 

= 5>>« Wa |A» a (B4) 

a 

= ^|«)Tr[4A] (B5) 



= (B6) 
The inverse of T a ^ w is given by 

= ^CT^W = (B7) 

T^J^)), = \A)) U . (B2) and hence T CT ^ W is unitary. 



